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nium conversion was established at P M j^ < 8.2 • 10 11 (90%C.L.) in 0.1 T 
magnetic field, which implies consequences for speculative extensions to the 
standard model. Coupling parameters in R-parity violating supersymmetry 
and the mass of a flavour diagonal bileptonic gauge boson can be significantly 
restricted. A Zs model with radiative mass generation through heavy lepton 
seed and the minimal version of 331-GUT models are ruled out. 
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A new upper limit for the probability of spontaneous muonium to antimuo- 
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At present all confirmed experimental experience is in agreement with conserved lepton 
numbers. Several solely empirical laws appear to hold simultaneously including multiplica- 
tive and additive schemes [1]. No associated symmetry has yet been identified thus leaving 
lepton numbers in a unique status in physics, since flavour mixing in the quark sector is well 
established and described by the Cabbibo-Kobayashi-Maskawa matrix. The standard model 
in particle physics assumes additive lepton family number conservation, and any observed 
violation would be a clear indication of new physics. In many speculative theories, which 
extend the standard model in order to explain some of its features such as parity violation 
in the weak interactions or CP violation, lepton flavours are not conserved. These theories 
have motivated a variety of dedicated sensitive searches for rare decay modes of muons and 
kaons [2] and for neutrino oscillations. 

Of particular interest is the muonium atom (M=/x + e~) which consists of two leptons 
from different generations. As the electromagnetic part of the binding is well described 
by electroweak standard theory it renders the possibility to determine accurate values of 
fundamental constants as well as searching for additional, yet unrevealed electron-muon 
interactions. A spontaneous conversion of muonium into antimuonium (M = /i~e + ) would 
violate additive lepton family number conservation by two units, however, it is allowed by a 
multiplicative law. This process could play a decisive role in many speculative models (Fig. 
1) [4,3,5-9]. 

The measurements reported here were performed with the Muonium - Antimuonium - 
Conversion - Spectrometer (MACS) whose design is based on the observation of M atoms 
in vacuo. In matter the possible conversion is strongly suppressed mainly due to the loss 
of symmetry between M and M due to the possibility of /i~ transfer in collisions involving 
M [10,11]. The required signature of a conversion process is the coincident identification of 
both the electron and positron released in the decay of the antiatom [12,13]. An energetic 
electron (e _ ) arises from the decay — > e~ + + V e with a characteristic Michel energy 
distribution extending to 53 MeV [14], and a positron (e + ) appears with an average kinetic 
energy of 13.5 eV corresponding to its momentum distribution in the atomic Is state of M 
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[15]. 

The setup has a large acceptance for these charged final state particles (Fig. 2). Its 
symmetry for detecting M and M decays through reversing all electric and magnetic fields 
is exploited in regular measurements of the M atom production yield which is required 
for normalization and, in addition, for monitoring detector performance. As a particular 
advantage, systematic uncertainties arising from corrections for efficiencies and acceptances 
of various detector components cancel out. 

The experiment utilizes the world's brightest continuous surface muon channel 7rE5 [16] 
at the Paul Scherrer Institut (PSI) in Villigen, Switzerland. It provides a central momentum 
p=26 MeV/c, a momentum bite Ap/p=5% and rates up to 8- I0 e fi + /s. The beam passes 
through a 280 /im scintillation counter and a 270 /im mylar degrader. Muonium atoms are 
formed by electron capture with 60% yield after stopping the /i + in a Si0 2 powder target of 
thickness 8 mg/cm 2 and supported in vacuo by a 25 /im aluminum foil with 30° inclination 
with respect to the muon beam axis. Most of the atoms emerge from the powder grains 
into the intergranular voids and then on average 3.3 % of them leave the target surface with 
thermal Maxwell-Boltzmann velocity distribution at 300 K [17]. 

When searching for M decays the energetic e" is detected in a magnetic spectrometer 
operated at 0.1 T magnetic field and covering 0.73*47r solid angle around the M production 
target. It has five concentric multiwire proportional chambers with radii of 8.2 to 32.0 cm 
and active lengths of 38 to 80 cm. They all are equipped with two planes of segmented 
helical cathode stripes for measuring radial, angular and axial coordinates. The efficiency 
for full particle track reconstruction is 82(2) % and 90(2)% if axial information is not re- 
quested. The momentum resolution at 50 MeV/c is 54(2) % yielding a probability of 10~ 5 
for misidentifying the charge of the particle. It is limited by the 2 mm wire spacing in the 
proportional chambers. The chambers are surrounded by a 64-fold segmented hodoscope 
providing a time resolution of better than 1 ns. Subsequent to the \i~ decay the atomic shell 
e + is accelerated to typically 7 keV in a two stage electrostatic device. It is guided in an 
axial 0.1 T magnetic field along a 3 m long transport region onto a microchannel plate de- 
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tector (MCP) with resistive anode readout and 0.55(3) mm intrinsic position resolution. A 
35 mg/cm 2 magnesium oxide coated carbon foil in front of a commercial MCP triplet stack 
serves as an additional secondary electron emitter and yields a fourfold enhancement of the 
detection efficiency to 64(2)% and reduces background counts of low energy ions trapped in 
the magnetic field [18]. The transport system is momentum selective due to a 90° horizontal 
bend of radius 35 cm and a collimator consisting of 40 cm long, 1 mm thick and 9 mm sep- 
arated copper sheats which act to suppress particles with longitudinal momenta exceeding 
750 keV/c because their gyration radii exceed 4.5 mm in the magnetic guiding field. The 
field gradient in the bend region causes a vertical drift for charged particles proportional to 
their momenta. It is compensated for 7 keV e + by a 270 V/cm transverse electrostatic field 
region of length 40 cm preceding the bend which also deflects low energy ji + and ions. 

Positrons are uniquely identified by annihilation radiation when striking the MCP which 
is centered inside of a barrel shaped 12-fold segmented barrel shaped pure Csl crystal de- 
tector. Using e + from radiactive sources the acceptance for at least one of two 511 keV 
annihilation photons was determined to be 79(4)% with 4.5(3)ns (FWHM) time resolution 
and 350(20) keV (FWHM) energy resolution. The transport path has 80(2)% transmission 
and conserves transverse spatial information of the decay vertex. It can be reconstructed 
radially to 8.0(4) mm and axially to 8.6(5) mm (FWHM), if in addition track parameters are 
used from the energetic e~ detected in the wire chambers. The primary limitations on the 
resolution arise at high energies from proportional chamber wire spacing and at low energies 
from multiple scattering in the 1 mm carbon fiber beam tube. 

During data taking, every 5 hours the M production yield was determined at low beam 
rates (2-10 4 // + /s) using a method which is based on a model established in independent 
dedicated experiments [17]. It assumes M production inside of the target where the /z + 
stop. Their escape into vacuum is governed by a random walk. The number of atoms in 
the fiducial volume was determined mainly from the distribution of time intervals between a 
beam counter signal from the incoming muon and the detection of the atomic electron on the 
MCP (Fig. 3a). The exponential background to these spectra was investigated using an Al 
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foil target and also negative muons stopped in Si0 2 powder (Fig. 3b). In both cases there is 
no M production but secondary electrons are released in a similar fashion as when a /i + decays 
in the target without M formation. The procedure is in agreement with results obtained 
from a maximum likelihood fit of the time distributions of M decay positions in four separate 
spatial regions above the production target. Due to finite acceptances and efficiencies of the 
detector components 5.0(2) • 1CT 3 of the incoming /z + on average were observed to decay 
as M atoms. The M atom fraction in vacuum per incoming /i + yields 2.0(1)% which is a 
factor of 2.5 below the results of earlier dedicated experiments where higher fractions were 
achieved due to narrower beam momentum bites [17]. The SiC>2 targets were replaced twice 
a week, since the M yield deteriorated on a time scale of a week associated with the release 
of H 2 molecules from the powder. 

The final search result was obtained in three data taking periods with a total duration 
of 6 month distributed over 4 years (Table I) during which the sensitivity of the instrument 
was constantly improved. In total 5.6(1)-10 10 M atoms were investigated for M decays. 
Two major sources of potential background were identified: (i) accidental coincidences of 
energetic e~ produced by Bhabha scattering of e + from M decays and scattered e + on the 
MCP and (ii) the allowed muon decay mode /j, —> e + e + e~z/ e 77 M with branching ratio 3.4- 10~ 5 
which could release an energetic e~ and a low energy e + which are dtected and the second e + 
escapes unobserved. As the expected e + flight time was determined in M decays to 76(1) ns 
with a 3.3(1) ns wide (FWHM) distribution, a narrow coincidence window of ±4.5 ns was 
applied to suppress both processes. Momenta of energetic e~ were accepted between 15 
and 90 MeV/c, which cuts 50(3)% of the Bhabha scattering events but only 8(2)% of the 
Michel spectrum contents. For the rare muon decay the bend of the e + transport system 
spoils the vertex reconstruction because of the significantly higher e + momenta. Positrons 
were required to deposit between 0.2 and 1.0 MeV in the Csl calorimeter and to be detected 
within ±6 ns of a hit on the MCP. The decay vertex had to be reconstructed within ±12 mm 
radially and was required to be not more than 3 cm upstream of the SiC>2 target. 

There was one event which passed all these required criteria within 3 standard devia- 

5 



tions of each of the relevant distributions which were established in the regular M control 
measurements prior to processing the search data (Fig. 4). The expected background due 
to accidental coincidences is 1.7(2) events and was evaluated using data with a TOF 4.5 ns 
above the expected value. 

The finite fiducial volume of diameter 9 cm and length 10 cm and the finite observation 
time are taken into account by a correction factor of 0.82(2), which was determined in a 
Monte Carlo simulation of the conversion process. As a result an upper limit on the con- 
version probability of P M m < 8-2 • 10 _11 /Sb (90% C.L.) was found. The factor Sb describes 
the suppression of the conversion in external magnetic fields due to the removal of degen- 
eracy between corresponding levels in M and M. It depends on the interaction type (Table 
II). The reduction is strongest for (V±A)x(V±A) [20,19]. For these cases the traditionally 
quoted upper limit for the coupling constant is 



where Gp is the weak interaction Fermi coupling constant. 

This new result which exceeds previous bounds [13] by a factor of 35 has some impact 
on speculative models. A proposed Z 8 model is ruled out with more than 4 generations of 
particles where masses could be generated radiatively with heavy lepton seeding [5]. 

A new lower limit of m x ±± > 2.6 TeV/c 2 *g 3 i (95% C.L.) on the masses of flavour diagonal 
bileptonic gauge bosons is found in GUT models, which is well beyond the value extracted 
from direct searches, measurements of the muon magnetic anomaly or high energy Bhabha 
scattering [8,21], with g 3t of order unity and depending on the details of the underlying 
symmetry. For 331 models this translates into mx±± > 850 GeV/c 2 , which excludes their 
minimal version in which an upper bound of 800 GeV/c 2 has been extracted from electroweak 
parameters [9,22]. The 331 model is only viable in an extended form with, e.g. a Higgs octet 



In the framework of R-parity violating supersymmetry [6,7] the bound on the coupling 




< 3.0 • HT 3 G F (90%aL.) 



[23] 
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parameters could be lowered by a factor of 15 to | A^A^ |< 3 * 1CT 4 for assumed super- 
partner masses of order 100 GeV/c 2 [24]. Further the achieved level of sensitivity allows 
to narrow slightly the interval of allowed heavy muon neutrino masses in minimal left-right 
symmetry [4] (where a lower bound on G M ^ exists, if muon neutrinos are heavier than 35 
keV/c 2 ) to 40 keV/c 2 up to the present experimental bound of 170 keV/c 2 [25]. 

In minimal left right symmetric models, in which MM conversion is allowed, the process 
is intimately connected to the lepton family number violating muon decay /i + — > e + '+v fl +u e . 
Within these models the limit achieved in this experiment this decay is not an option for 
explaining the excess neutrino counts in the LSND neutrino experiment at Los Alamos [26]. 

The consequences for atomic physics of M are such that the expected level splitting in 
the ground state due to M — M interaction is below 1.5 }^z/^/S^ reassuring the validity of 
recent determinations of fundamental constants from atomic spectroscopy of the atom. 

Further increased sensitivity to MM conversion could be expected, if the M production 
efficiency could be enhanced, or from muon sources with significantly higher fluxes, e.g. at 
the front end of a future muon collider. 
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TABLES 

TABLE I. Summary of three data taking series. The number of produced muonium atoms and 
the upper limits on the conversion probability in 0.1 T magnetic field are given. The probability 
in zero field can be obtained through division by the interaction dependent reduction factor Sb of 
table II. 



measurement duration of 


number of 


P mm(0-1 T) 


data taking 


muonium atoms 


(90%C.L.) 


1 210 h 


1.4(1) • 10 9 


2.8 • 10~ 9 


2 230 h 


2.3(2) • 10 9 


1.8 • 10~ 9 


3 1290 h 


5.2(1) • 10 10 


8.9 • HT 11 


all 1600 h 


5.6(1) • 10 10 


8.2 • 10- 11 


TABLE II. Magnetic field correction factor Sb of muonium- antimuonium conversion probabil- 


ity for muonium atoms with statistically populated ground states [19,20]. 




interaction type 


2.6 fxT 0.1 T 


100 T 


S S 


0.75 0.50 


0.50 


P P 


1.0 0.9 


0.50 


(V±A)x(V±A) or (S±P)x(S±P) 


0.75 0.35 


0.0 


(V±A)x(VtA) or (S±P)x(StP) 


0.95 0.78 


0.67 
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FIGURES 



mmb_f 1 .ps 

FIG. 1. Muonium-antimuonium conversion in theories beyond the standard model. The 
interaction could be mediated, e.g. by (a) doubly charged Higgs boson A ++ [3,4], (b) heavy 
Majorana neutrinos [3], (c) a neutral scalar <E>tv [5], e.g. a supersymmetric r-sneutrino v T [6,7], or 
(d) a bileptonic flavour diagonal gauge boson X ++ [8,9]. 




FIG. 2. The MACS apparatus at PSI searching for muonium-antimuonium conversion. The 
signature requires observation of the energetic e~ from /i~ decay of M in a magnetic spectrometer, 
in coincidence with the atomic shell e + . The latter is accelerated and magnetically guided onto 
a microchannel plate (MCP). In addition at least one annihilation photon is required in a Csl 
calorimeter. 
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mmb_f 3 .ps 



FIG. 3. Muonium production was continually monitored using the characteristic time of 
flight distribution of atomic electrons (a). The indicated exponential background was verified by 
demonstrating that there is only an exponentially decaying background for negative muons on SiC>2 
powder (b). 



mmb_f4a.ps 



mmb_f 4b .ps 



FIG. 4. The distribution of the distance of closest approach (Rrf ca ) between a track from 
an energetic particle in the magnetic spectrometer and the back projection of the position on 
the MCP detector versus the time of flight (TOF) of the atomic shell particle, for a muonium 
measurement (left) and for all data recorded in the third data taking period of 1290 h while 
searching for antimuonium (right). One single event falls within a 3 standard deviations region of 
the expected TOF and Rdca indicated by the ellipse. The events concentrated at low TOF and low 
R^ ca correspond to the allowed decay \i — > 3e + 2u. 
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